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ABSTRACT 

We report the first detection of thermal water vapor emission in the 557 
GHz, 1 iq — loi ground state transition of ortho-^O toward VY Canis Majoris. In 
observations obtained with the Submillimeter Wave Astronomy Satellite (SWAS), 
we measured a flux of ~ 450 Jy, in a spectrally resolved line centered on a velocity 
v lsr — 25 km s' 1 with a full width half maximum of ~ 35 kms -1 , somewhat 
dependent on the assumed line shape. We analyze the line shape in the context of 
three different radial outflow models for which we provide analytical expressions. 
We also detected a weaker 557 GHz emission line from W Hydrae. We find that 
these and other H 2 0 emission line strengths scale as suggested by Zubko and 
Elitzur (2000). 

Subject headings: circumstellar matter - infrared: stars - stars: mass loss - 
stars: winds, outflows- stars: individual (VY Canis Majoris, W Hydrae) 


1. Introduction 

VY Canis Majoris is a highly luminous, variable supergiant. It is strongly obscured; 
only 1 % of the total luminosity is detected at optical wavelengths. The line of sight 
optical depths to the star at wavelengths of 1.65, 2.26 and 3.08 /mi are 1.86 ± 0.42, 0.85 
± 0.20 and 0.44 ± 0.11 (Monnier et al. 1999). Correspondingly, the star exhibits a huge 
infrared excess. The star’s distance is believed to be D = 1500 pc, implying a luminosity 
L = 5 x 1O 5 L 0 . (cf. Lada & Reid 1978, Richards et al. 1998). A high signal-to-noise ratio 
infrared spectrum shows the star enveloped in an optically thick dust cloud whose inner 
edge is at a temperature of several hundred degrees (Harwit et al. 2001, hereafter Paper 
I). At greater distances from the star the spectrum is dominated by amorphous silicate 
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features. 

Strong outflows have been reported by Buhl et al. (1975) who detected thermal 
emission in the vibrational ground state (v = 0, J = 2 - 1) 86.84686 GHz SiO line. Reid and 
Dickinson (1976) interpreted these in terms of a stellar radial velocity of 17.6 ± 1.5kms _1 
with respect to the local standard of rest, and an expansion velocity of 36.7±2.0 kmsW 
Three thermally emitted mid-infrared water vapor emission lines seen in high-resolution 
spectra obtained with the Short Wavelength Spectrometer (SWS) on the Infrared Space 
Observatory (ISO) indicate a mean radial velocity of order 20 ± 2kms _1 and a 25 km s^ 1 
outflow velocity (Neufeld et al., 1999). These water vapor outflow velocities are significantly 
lower than the velocities for SiO cited by Reid & Dickinson (1976) or the 32kms _1 velocities 
that Reid & Muhlcman (1978) reported for 1612 MHz OH maser outflow. 

The effective temperature of the central star appears to be T* = 2, 800 K (Monnicr et 
al. 1999). The combination of effective temperature, distance, and inferred luminosity fixes 
the star’s radius at R* ~ 15 AU. In the near-infrared at 0.8, 1.28 and 2.17/mi, VY CMa 
appears embedded in an elongated circumstellar envelope with respective dimensions of 
60 x 83, 80 x 166 and 138 x 205 mas (Wittkowski et al. 1998). Working at 11 /mi, Danchi et 
al. (1994) reported a variable photospheric radius ranging from 9.5 to 11 mas, where 10 mas 
corresponds to a stellar radius of 15 AU; the surrounding dust shell had a corresponding 
inner radius of 40 to 50 mas and an outer radius of 2.5 arcsec. 

To model the dust disk surrounding VY CMa, Efstathiou & Rowan-Robinson (1990, 
hereafter ERR) assumed that the central star is orbited by a massive dusty disk, which 
flares out from the star at angles ±45° above and below the disk’s central plane. At angles 
> 45° off the plane, radiation escapes into space. ERR believe that the line of sight from 
Earth is inclined at about 43° to the disk’s orbital plane, so that we are looking almost 
tangentially through the outer layers of the disk. This radiative transfer model has been 
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widely accepted because it accounts for the dust spectrum, the high reddening, and the lack 
of visible radiation from the star. However, as Ivezic & Elitzur (1995) have emphasized, 
a correct model must self-consistently couple the equations of radiative transfer to the 
outflow dynamics. Zubko & Elitzur (2000, hereafter ZE) have recently carried out such a 
calculation, and used it to extensively model the semi-regular M7.5 red giant W Hydrae. 

In this letter we re-examine the outflow problem. In §2, we present new data on the 
strength and shape of the 557 GHz water vapor line observed with the Submillimeter Wave 
Astronomy Satellite (SWAS) in VY CMa and W Hydrae. In §3 we introduce complementing 
ISO data. §4 attempts to model the geometry of the VY CMa outflow, while §5 portrays 
its radial dependence. A concluding section summarizes our findings. 


2. Observations 

In the course of many hundreds of SWAS satellite orbits, between November 1999 and 
May 2000, we carried out observations on VY CMa in the 556.936 GHz 1 10 — l 0 i ground 
state transition of ortho water vapor. All observations reported here were made with the 
spacecraft in the nodding mode (Melnick et al. 2000), with the telescope’s primary held of 
view centered on a(2000) = 07 h 22 m 58.3 s , <5(2000) = —25° 46' 03'.'0. The observing time 
accumulated on source was 7902 minutes (131.7 hours) with an equal time devoted to a 
blank held displaced south in declination by 0.24°. 

The observed water vapor line, shown as a histogram in Figure 1, is centered at 556.89 
GHz, corresponding to a velocity with respect to the local standard of rest vlsr = 25kms _1 . 
A ht to the line shape discussed in §4 yields a full width at half maximum (FWHM) of 
~ 35 km s' 1 ; the peak signal to noise ratio is S/N — 7. The area under the ht divided by 


its FWHM yields an antenna temperature T\ = 0.033 ± 0.005 K ~ 450 Jy. 
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Using a beam centered on a(2000) = 13 h 49 m 02.1 s , 5(2000) = —28° 22' 03'.'0, we 
carried out observations on W Hydrae in January and February 1999, July 1999, July 2000, 
and February 2001 for a total of 126.7 hours, and also detected the 557 GHz line. The 
antenna temperature was T* A = 0.016 ± 0.004 K with a line width ~ 10.4 ±2.1 kms -1 . 

3. Complementary Water Vapor Data 

Neufcld et al. (1999) reported line strengths for more than forty thermal water vapor 
emission lines from VY CMa observed with ISO, in the wavelength range of ~ 30 to 45 /mi 
Three transitions, respectively at 29.84, 31.77 and 40.69 /an, were also observed with SWS 
in its Fabry-Perot mode and exhibited P Cygni profiles and emission lines peaked at an 
average velocity of vlsr ~ 20kms -1 , consistent with the SWAS observations obtained at 
a higher velocity resolution of ~ 1km s -1 . The outflow velocity of ~ 25kms -1 inferred by 
Neufcld et al. (1999), is also in reasonable agreement with our 557 GHz results. The line 
profile obtained with SWAS does not exhibit P Cygni profiles, nor would we expect it to 
because the continuum is too weak to be detected. Archival ISO data obtained with the 
Long Wavelength Spectrometer, LWS, provide additional flux levels for transitions at 66.44, 
108.078, 174.62, and 179.53 /mi. Table 1 compares these data on VY CMa to measurements 
on W Hydrae compiled by ZE and found to be in reasonably good agreement with their 
theoretical model for the star. The table shows the flux levels in VY CMa and W Hydrae 
to be approximately in the ratio 10:1. The VY CMa data scaled down by a factor of 10 
correspond well to the ZE model. 
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4. Modeling the VY CMa Outflow 


The outflow from red supergiants is generally assumed to be propelled by radiatively 
accelerated dust grains that transfer their momentum to ambient gas. The VY CMa mass 
loss rates cited by different observers range from 1 to 3x10 4 M q yr 1 (cf. Paper I for 
a summary). The water vapor outflow can be expected to be roughly 1% of this by mass 
and, given the high Einstein coefficient for the 557 GHz transition, A = 3.45 x 10 -3 s -1 , 
the observed emission has to be optically thick. Velocity gradients, however, permit the 
radiation to escape; otherwise we would see only a velocity component from the approaching 
part of the outflow, rather than from its entirety. 

Despite the disk-shaped outflow postulated by ERR, the structure of the cloud 
surrounding VY CMa is quite uncertain. Planetary nebulae often light up bipolar flows that 
may have been ejected during a previous evolutionary phase. Fig 2(a) provides coordinates 
for an axially symmetric outflow whose axis of symmetry is inclined toward the observer at 
an angle 9 0 . A spherically symmetric, constant velocity, optically-thick outflow (model A) 
would yield a parabolic line shape shown as the solid curve in Figure 2b. 


For the bi-conical geometry shown in Figure 2a, the flow directed at angle 6 to the 
observer’s line of sight is radially outward from the origin and describes a circular arc of 
radius aa! = aa". Its length l spans points a! and a" on a circle defined by the angle 9 0 , 
by the opening angle a of the cone, and by the tilt angle toward the observer, 9, of the 
conical surface joining the origin to points on the circular arc. The length of the arc, t can 
be expressed as 


= R sin 9 


7r — 2 sm 


sm ti 


( 1 ) 


The segment formed by joining the origin to points along this arc presents to the 
observer a pie-shaped area projected on the plane of the sky, of measure (R£ sind)/2 defined 
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by points a, a', and a". The projected velocity along the observer’s line of sight throughout 
this surface is v = Uocos 9; for an optically thick line the velocity distribution yielding the 
line shape is given by the relative values of (RP smd)/2 at different angles 9 = cos^ 1 (u/uo)- 

An optically-thick, bi-conical, constant-velocity outflow, symmetric about an axis 
perpendicular to the observer’s line of sight (i.e. with 9 0 = vr/2) and limited to angles less 
than a = 7t/3 from the axis of symmetry — with a density independent of angular distance 
from the axis (model B), yields a line profile shown by the dashed curve in Figure 2b. A 
uniform bi-conical outflow within an angle of a = 7t/4 from the axis, and with 9q = 47° 
(model C), would yield the bi-lobed dotted line shape shown in Figure 2b. Neither of these 
two models produces the line shape observed. 

To obtain the line profile for an optically-thick disk, we subtract these respective 
bi-conical outflows from a spherical outflow. For a disk with axis perpendicular to the 
line of sight to Earth, and a latitude-independent outflow into a disk with half angle ±30° 
above and below the symmetry plane, we can subtract the dashed curve in Fig 2b from the 
solid curve — i.e. models (A - B) — to obtain the fit shown by the dashed line in Fig 1. 
The line shape expected from the model proposed by ERR = models (A — C), would result 
from a subtraction of the dotted line in Fig. 2b from the solid curve. Figure 1 shows that 
models A and (A-B), normalized to the line strength and width, give rather better fits to 
the observed 557 GHz profile than does the ERR model (A-C). 


5. The Temperature/Radius Profile 

Zubko and Elitzur (2000) have devised the most comprehensive computational model 
of outflow from highly evolved stars developed to date. It derives a profile for temperature 
T and radial distance R from the star, T 2 R ~ constant. The gas in the outflow is primarily 
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heated through collisions with grains, cooled by H 2 O emission and adiabatic expansion, 
and effectively acts as though it was in thermal balance with the parent star throughout 
the flow. In apparent agreement with the general scaling properties proposed by Ivezic & 
Elitzur (1997) as applied to dusty gaseous outflows, Table 1 shows the good fit of the ZE 
model to the observed H 2 O line strengths in W Hydrae, that also scales well to VY CMa. 


Given this observed agreement, we adopt the ZE profile, T 2 R = constant, for VY 
CMa. Paper I suggests the existence of a dust photosphere at a temperature T p ~ 400 K 
in thermal equilibrium with the star. R p is a characteristic radius for this photosphere in 
thermal balance with the star, R 2 T p = R 2 Tf, so that the ZE profile can be taken to extend 
all the way to the stellar surface: T 2 R = T 2 R *. In Table 3, we list the radial distances at 
which the VY CMa water vapor lines emanate, on the assumption that the emitting region 
is optically thick - an assumption justified by the deep P Cyni profiles of even the highest 
rotational lines observed. The flux received from a circumstellar cloud subtending a solid 
angle ff at the observer is then 


F(u) 


kT o 2 kT R 
~ ~ ^D 2 


2u 2 


( 2 ) 


where we assume that all the radiation is in the Rayleigh Jeans part of the spectrum. Table 
3 shows this assumption to be consistent; D is the distance to the star, R is the radial 
distance of the emitting H 2 0 region from the star, and ^ is a constant of order unity that 
depends on the geometry of the outflow. For a sphere it has a higher value than for a 
flattened disk. Table 3 assumes ifj ~ 1. Substituting for the stellar parameters T* = 2800 K, 
R* = 2.25 x 10 14 cm, D = 4.5 x 10 21 cm and taking the emission to be blackbody in the 
frequency interval observed, we can solve for the characteristic distances from the central 
star and the corresponding temperatures at which the different spectral lines appear to be 
emitted. 


R 


1.75 x 10 25 F(z/)A 2 i 2 /3 




( 3 ) 
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where F(y) is given in Janskys, and where R and A both are in units of centimeters. The 
calculated distances for the last three lines in Table 3 use the full emission line strengths 
corrected for P Cygni absorption of the continuum. For the longer wavelength lines observed 
with the ISO LWS, the resolution is insufficient to determine whether self-absorption 
weakened the lines. 


6. Conclusions 

We have obtained a high-resolution spectrum of the 557 GHz water vapor emission 
line from the circumstcllar cloud surrounding VY CMa, and have obtained a weaker signal 
from W Hydrae. The observed flux for VY CMa is ~ 450 Jy. The line shape is symmetric 
around a recession velocity vlsr — 25kms~ 1 and has a FWHM of order ~ 35kms~ 1 . We 
have developed analytical line fits for several plausible models of the outflow. Observations 
that will be possible with the Herschel space telescope should yield improved line shapes 
that will further constrain outflow models. Archival and published ISO data for VY CMa 
and W Hydrae, suggest the adoption of a model developed by Zubko & Elitzur (2000) to 
determine the distance from the central star and the temperature regime around VY CMa 
from which the H 2 0 lines emanate. 
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Figure Captions 

Figure 1. The 557 GHz water vapor spectrum of VY CMa obtained with SWAS. The histogram 
gives the SWAS observations. The solid line is a fit for a spherically symmetric outflow 
(model A). The dashed line is for a disk viewed edge-on (models A - B); the dotted curve is 
a fit for the ERR outflow (models A - C). All three curves are normalized to the observed 
line profile (see text and Figure 2 for details). 

Figure 2. (a) The geometry of the outflow for which the velocity profiles were derived. R is the 
radius of a sphere whose radial extent is that of the outflow. For a bi-conical or disk-shaped 
flow, the angle 6o gives the tilt angle toward the line of sight. Throughout, we assume for 
simplicity that the symmetry axis and the line of sight define a plane of mirror symmetry 
in the outflow. The angle a is the opening angle of the conically shaped flow as measured 
from the symmetry axis. The angle 6 defines the radial direction of motion of a volume in 
the flow relative to the line of sight to the observer. Points a, a', and a" lie in the plane of 
the sky. (b) The solid curve gives the expected line shape from a spherical outflow (model 
A). The dashed curve gives the expected line shape for a bi-conical outflow (model B) 
with 0 O = 90° and a = 60°, while the dotted line (model C) gives the line shape for an 
inclined bi-conical outflow with 6 0 = 47° and a = 45°. Subtracting the bi-conical outflows 
from that for spherical outflows, respectively, yields the dashed and dotted line prohles for 
disk-shaped outflows in Figure 1. (See text for details). 
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Table 1. Relative 0 -H 2 O Flux Levels from VY CMa and W Hydrae 


Wavelength 

A (/urn) 

Transition 

W Hydrae a 

10 _20 W cm -2 

Model 6 

10 -20 W cm -2 

VY CMa c 

10~ 19 W cm -2 

538.289 

lio —*■ loi 

0.45 

0.76 

0.34 

179.527 

2l2 —> loi 

8.66 

12.48 

8.7 

174.624 

3o3 —► 2 12 

9.21 

8.31 

00 

bo 

108.073 

2 l 2 —> lio 

13.00 

15.47 

13.2 

66.438 

330 —> 221 

22.90 

18.45 

17.2 

40.691 d 

432 —> 3o3 

23.00 

36.07 

36.9 

31.772 d 

441 —> 3l2 

63.0 

39.95 

24.1 

29.837 

725 —► 616 

32.00 

30.72 

23.5 


a ISO from Neufeld et al. (1996), Barlow et al. (1996); SWAS data, this paper 
b W Hydrae model from ZE, and 538 fim flux through 
private communication from Drs. Zubko Sz Elitzur 
c For sources of data see text 
d Blending may need to be taken into account 


Table 2. Characteristics of VY CMa and W Hydrae 


Star 

VY CMa a 

W Hydrae 6 

Star’s Distance (pc) 

1500 

115 

Star’s Luminosity (Lq) 

5 x 10 5 

11,050 

Star’s Temperature (K) 

2800 

2500 

Star’s Radius (cm) 

2.25 x 10 14 

3.9 x 10 13 

Final Outflow Velocity (km s —1 ) 

20 

10 

Optical Depth 

Aj ~ 3.2 

0.83 at 5500A 

Gas-to-dust ratio 

~ 100 

850 

Mass loss rate Mq yr -1 

~2x 10“ 4 

2.3 x 10“ 6 


a Data from Neufeld et al. (1996), Barlow et al. (1996) 
b For sources of data see text 




















Table 3. Ortho-E^O Emission: Characteristic Distances and Temperatures 


Wavelength 

(fnn) 

Flux 

(Jy) 

Radius 0 

(10 16 cm) 

Temperature* 2 

K 

rp b 

1 U 

K 

Ti b 

K 

538.289 

450 

8.1 

148 

60.95 

34.23 

179.527 

3490 

7.3 

155 

114.35 

34.23 

174.624 

3440 

7.2 

156 

196.72 

114.35 

108.073 

3190 

6.9 

160 

194.05 

60.95 

66.438 

2560 

1.6 

332 

410.55 

194.05 

40.691° 

3360 

0.99 

431 

550.30 

196.72 

31.772° 

1710 

0.45 

626 

702.21 

249.37 

29.837 

1570 

0.39 

672 

1125.55 

643.34 


a Derived from constancy of TjR 1 / 2 , see text 
b Level energies divided by Boltzmann constant 
c Blending may need to be taken into account 
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